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1. Introduction 

The broad research programme of ultra-relativistic heavy-ion collisions aims at 
investigating a new state of deconfined strongly interacting matter. Experimentally, 
the hints for such a state, dubbed quark-gluon plasma, have been accumulated in 
investigations at CERN-SPS [Ij and further consohdated at BNL-RHIC [2]. While early 
ideas have been guided by asymptotic freedom considering the quark-gluon plasma as an 
ensemble of weakly interacting quarks and gluons, the paradigm has changed now to a 
strongly coupled quark-gluon plasma [3] , as enforced by the success of hydrodynamical 
concepts which point to an extremely fast thermalization [3] and an exceedingly low 
viscosity [H [5] . 

Parallel to the intense experimental efforts, which will soon proceed to a new era 
at CERN-LHC, much progress has been achieved on the theoretical side. In particular, 
first-principle (lattice) QCD calculations, directly based on numerical evaluations of the 
QCD partition function, are progressing and deliver information on the equation of state 
and related thermodynamic quantities such as susceptibilities etc. [6j. While final results 
on bulk properties of deconfined strongly interacting matter at finite temperatures 
seem to be achievable in the near future, the microscopic nature - with respect to 
the excitations in the considered medium - received less attention hitherto. 

In fact, a variety of phenomenological approaches basing on different microscopic 
pictures account fairly well for the bulk information obtained in lattice QCD. For 
instance, various effective quasi-particle models [3 El El [IDl [HI [121 [E] rely on 
qualitatively formulated / postulated excitations in the quark-gluon medium in order 
to achieve a parametrization of the equation of state adjusted to lattice QCD results. 
Certain symmetries of QCD are the basis of PNJL models (cf. |14|) which also need 
adjustments to lattice QCD data. Other approaches assume that the QGP may 
be described by a classical non-relativistic plasma of colour charges with Coulomb 
interaction [15] for temperatures Tc < T < ST^, where T^. is the pseudo-critical or 
deconfinement temperature. Analytical attempts to derive the equation of state from 
QCD by perturbative means [16] are in agreement with the non-perturbative lattice 
QCD results for temperatures T > ST^. These approaches also provide some legitimation 
for a quasi-particle picture that goes beyond strict perturbation theory [T7] . 

At asymptotically large temperatures suitable techniques like dimensional reduction 
|19] bring further insight into the structure of the theory and are applicable, e.g., to the 
physics of the early universe. (Here, at temperatures 1 - 100 GeV, the heavy quark 
sector becomes important [iH]-) Strict perturbative expansions have been pushed to 
order g^lng (see [T9l [20] and references therein) and slightly beyond [18], where g is 
the strong coupling. The successful hard thermal loop/hard dense loop (HTL/HDL) 
resummation scheme pi] employs HTL self-energies, which correspond to the high 
temperature/density limit of one- loop self-energies, thus, setting all quark masses to 
zero. The systematics of one-loop and HTL self-energies, as well as the corresponding 
dispersion relations were studied in [22] for zero quark masses. Further steps towards 
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going beyond the one-loop approximations have been attempted in [23], where a Dyson- 
Schwinger scheme is set up in ladder approximation and in the chiral limit focussing on 
the fermion spectral function. 

It is the temperature region — ST^ which is of utmost relevance for heavy-ion 
collisions. Here, the information on the excitation spectrum of deconfined matter is fairly 
scarce. In [2l], the poles of quark and gluon propagators in Coulomb gauge have been 
analyzed with the result that they may be parametrized by an energy [uj) - momentum 
{k) relation oj"^ = k'^+m'^ g with nig^g/T = g and ^q^g ranging from 1.2 till 3.9 depending 
on parton species and temperature. Karsch and Kitazawa [25] analyzed, in quenched 
approximation and in Landau gauge, spectral properties of quarks at 1.5Tc and 3Tc and 
at zero momentum as a function of the bare quark mass m. An important result of these 
investigations is the confirmation of a mass gap ruq = ^q{T)T with a weak temperature 
dependence in the function C,q{T). Another important result in [25] concerns the quark 
mass dependence of pole positions and residues of the quark propagator, which was 
found to be qualitatively different from the expected perturbative pattern. 

For the latter one, systematic investigations are hardly found in the literature. The 
thermal self-energies of gluons (e.g. [26]) and quarks (e.g. have been calculated for 
massive quarks. Fermion mass effects on the dispersion relations have been studied by 
several authors [271 [2H1 [291 SO] with the restriction to the long-wavelength limit of the 
dispersion relations or to small (soft) masses m < gT or m < gfi, where fi is the quark 
chemical potential. 

A more detailed investigation of the general mass dependence in the dispersion 
relations is desired, e.g., to uncover effects of masses m ~ 0(T), as employed in lattice 
QCD calculations [31]. A deeper understanding of the impact of heavy quark masses 
on the equation of state of strongly interacting matter relevant for the early universe is 
also of interest. 

The chiral extrapolation, and thus the quark mass dependence, is an important 
issue not only in thermo-field theory. Also in low energy effective field theories, such 
as chiral perturbation theory, the chiral extrapolation of quantities like the nucleon 
mass is of great interest [32]. In hot QCD it is particularly challenging as lattice QCD 
evaluations have been performed often for nonphysically heavy quarks due to technical 
limitations. Apart from the necessary continuum extrapolation, the extrapolation to 
physical quark masses is an inevitable step towards obtaining useable results. 

Given this motivation, we present here the systematics of the quark mass 
dependence of quasi-particle dispersion relations in hot one-loop QCD. We mention as 
further motivation that in a series of papers [H El [10] the successful description of lattice 
QCD results on bulk properties of strongly interacting matter was demonstrated within 
a phenomenological quasi-particle model employing approximate one-loop dispersion 
relations for quarks and gluons. The model [3, [9l [101 [H] goss beyond perturbation 
theory by using an effective coupling rather than the strong coupling. This motivates 
us not to restrict ourselves to the weak coupling regime. 

Our paper is organized as follows. In section [21 we present numerical results of 
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the thermal part of gluon one-loop self-energies and related quantities, in particular 
the dispersion relations which offer a glimpse on the relevant excitations. Section [3] 
illuminates the same quantities but now for quarks. Section H] is devoted to the gauge 
dependence of our results. For this purpose, we contrast our results with calculations 
in Coulomb gauge. Of course, there are regions in parameter space where our results 
coincide with the HTL approximation, which has been proved to be gauge independent 
in [33] . The summary can be found in section [5l Appendix A lists necessary relations 
for the chiral expansion of the asymptotic thermal masses and Appendix B contains 
the decomposition of the quark propagator with non-zero quark mass into different 
physical excitations according to the spinor structure of the propagator. In Appendix 
C, figures for quantifying the self-energies are collected. In Appendix D, we present 
an asymptotic Dyson-Schwinger approach to estimate higher order contributions to the 
mass dependence of the asymptotic thermal masses in the case of Abelian gauge theory. 

2. One-loop gluon excitations 

In the following, we consider the quark masses as well as the coupling g as external 
parameters. Our results are presented for two-flavor QCD {Nf = 2) with one 
independent quark mass parameter to limit the dimensions of parameter space. However, 
a generalization to cover a hierarchy of different quark masses is straightforward. We 
focus on the thermal parts of the self-energies at zero chemical potential (quark-anti- 
quark symmetric hot QCD medium), having in mind, however, the remarks given in 
[28]. 

2.1. Gluon self-energy 

The gluon one-loop self-energy may be decomposed as 



with Lorentz indices /i, u and colour indices a, b of the adjoint representation of the 
gauge group. The contributions of gluon loop {3g), gluon tadpole (4(7), and ghost loop 
{ghost) are grouped in the Yang-Mills contribution "YM" . The quark loops are labelled 
by "g" and the sum runs over all quark flavors included. They read as functions of the 
four-momentum K of the considered gluon 
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In the above expressions the fermion (F) and boson (B) propagators are the bare 
propagators S = So{K) = {m — fC)AF{K) and = {Vq)^^ = d^y^B^K) in Feynman 
gauge suppressing the exphcit colour indices for notational convenience, with the thermal 
propagators Af{K) = (cj^ + A;^ + m^)~^ and Ab{K) = (a;^ + with Matsubara 

frequencies Un = 2miT and uon = (2n + l)7rT for bosonic and fermionic degrees of 
freedom, respectively. The trace "tr" has to be taken over spinor and colour indices 
and P and Q = P — K are the internal loop four-momenta. For the expressions of the 
various vertices F we refer to standard textbooks on thermo-field theory, e.g. [31] . 
The tensor 11^^, consists of two independent scalar functions according to 

UL = C.n?. + ^J^ni, (4) 



where the utilized four-transverse projectors {K'^^J^i = 0) are given by 



^2 ' ' 
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with N^, = (K^iKu) - u^K^) {{Kuf - K^) and K^N^ = 0. The superscripts T 
and L indicate that the tensors and ^^^^ project on subspaces transverse and 
longitudinal to the three-momentum k, respectively. After performing the summation 
over the Matsubara frequencies, integrating over the angular part of d^p and analytically 
continuing into Minkowski space-time one obtains for the individual scalar functions [26] 
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Figure 1. The scaled real part of the quark loop contribution to the transverse gluon 
self-energy according to © as a function of the scaled energy for Nj = 2 degenerate 
quark flavors with mass m/T = (blue short-dashed curve), 0.5 (red solid curve) and 
0.7 (green long-dashed curve). The kinky structure aiuj — y/k'^ + Am? is clearly visible. 
The black dotted line represents ll? ~ k"^ . Thus, the intersection with the self-energy 
function would represent the sole solution of the dispersion relation if there were no 
other contributions. The curves are for g = 1 and k — gT. 



and ep = \Jp^ + m?. K ~ (to", k) consists of the two components uj = K ■ u and k 
with |k| = /c = A/(^^"^iip"^^~^ in the thermal medium with four-velocity u^. np^B 
are the standard Fermi and Bose distribution functions and the gauge group factors are 
Ca = N^, C2 = 1/2 for a SV{N^) gauge group with A^, = 3 for QCD. 

For illustration, figure [T] exhibits the real part of 11^ as a function of u for 
three different values of the quark mass at fixed k. The real parts of the quark 
loop contributions H'^ ^ develop a non-analytical (continuous but not continuously 
differentiable) behaviour at the threshold u"^ = k'^ + 4m^ for the production of a real 
quark-anti-quark pair. This behaviour is, however, related to the imaginary part of 
n|,, which is also non-zero above the threshold u"^ > k"^ + Am?. The imaginary parts 
in such one-loop calculations are known to have no direct physical meaning since the 
quasi-particle damping rates (related to the widths of quasi-particle peaks in the spectral 
function) turn out to be gauge dependent, both, in magnitude and even in sign. Braaten 
and Pisarski [21] pointed out the necessity of resumming diagrams of higher-loop order 
(HTL resummation) to get a definite result. However, the excitation energies of the 
quasi-particles considered in this paper are rather independent of the imaginary parts. 

With the replacements g e, m —>■ mg, Ca — > 0, C2 — > 1 and Nf 1 the 
quantities 11^ ^ would represent the thermal transverse and longitudinal one- loop photon 
self-energy contributions in a hot Abelian (QED) plasma, respectively. Such plasmas 
are presently of interest [35l|36], as their production under laboratory conditions seems 
feasible with high-intensity lasers. The temperature of such a plasma is expected to be in 
the range of T ~ 1 ... 10 MeV while the electron mass me ~ 511 keV sets an additional 
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energy scale. Thus, in this temperature region, in particular for lower temperatures, the 
ratio vrif,jT is not small, which might lead to modifications in certain plasma properties 
as compared to the ultra-relativistic case. Consequently, a non- negligible electron mass 
is expected to infiuence physical observables such as production rates of pions and muons 
[35] and of photons [37] . 

2.2. Gluon dispersion relations 

The dressed gluon propagator may be decomposed as 

V;l{K) = S^' (^^l^Ar + + P^) (12) 

with gauge fixing parameter p = 1 in Feynman gauge and the projectors ^^i/" are defined 
in the explanation below equation (jl]). The transverse and longitudinal propagator parts 
are related to the corresponding self-energies by Dyson's equation via 

" K^-Wr{uj,ky " k^-iiiiu^k) ^^^^ 

with K,l = ^t!1 + Y.,K,l- 

The gluon one-loop self-energies possess in general real and imaginary contributions. 
Here, we are interested in the dispersion relations of gluonic quasi-particle excitations 
with negligible damping rates, which are characterized by 

Re (A^i) = 0. (14) 

The gluon excitation energies, which are the real solutions of (|T4l) . may be represented 
by uj^^^ = + gT,Lml with = \g'^T'^{Nc + Nf/2). The contribution m^^T.L 
to the quasi-particle excitation energy encodes the effects of the thermal medium 
and represents a thermal mass. Therefore, the light cone is given by the planes 
Qt,l = in this representation. For the transverse gluon excitation one can write 
Qt = Re Uxiuj^k), k,m)/m'^ and for longitudinal gluons (plasmons) one finds Ql = 
K^Re ni(co'(fc), fc, m)/(fc^mp. In the presently considered case, the parameter space is 
three-dimensional for fixed temperature: The functions Qt,l depend on m, g and k. 
Note that the g dependence cannot be scaled out, unlike in the HTL approximation. 
We are interested in the regions m ~ T and ~ T. The first range is determined 
by the lattice calculations [31] which employ "lattice quark masses" m oc T; the 
interesting range is extended till the chiral limit m = 0. The second range is determined 
by the observation [7| that thermal excitations with A; ~ T essentially contribute to 
thermodynamic quantities. A survey of these functions is exhibited in figure [2] for 
Nf = 2 quarks of mass m. 

The overall observation is, that the thermal contribution to the transverse gluon 
mode stays above the light cone in the whole parameter range, whereas it drops 
and rapidly approaches the light cone for the longitudinal mode at sufficiently large 
momenta, independent of g and m. Both, and Ql, are monotonically decreasing 
with increasing mass m, indicating that heavy quarks, compared to massless quarks. 
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contribute not as much to the thermal mass of the gluon. The m dependence is, 
though, rather weak because the thermal masses are dominated by the strong gluonic 
self-interaction. For the transverse mode, the function Qt approaches unity in the limit 
— >• cx) and m — 0, exposing that nig is the asymptotic mass, indeed. We emphasize 
the strong g dependence when considering the region (7 ~ 1, see figure El (Note that 
g = 1 corresponds to as ~ 0.08, whereas g = 0.3 translates to as ~ 1/137 and g = 3 
means as ~ 0.72). Clearly, for large values oi g, higher order contributions are expected 
to become important. 

Such 3D plots are useful for surveys, however, cuts allow for a better quantitative 
representation of results. Corresponding figures are relegated to Appendix C (cf. figures 
[CT]and [C2|l . 

2.3. Gluon plasma frequency and asymptotic behaviour 

It is instructive to study some special limiting cases of the dispersion relations, where 
the mass dependence of the thermal mass part of the excitation energies can be given 
explicitly by analytical formulae. The gluon plasma frequency is defined as the long 
wavelength limit of the gluon dispersion relation Upi = limu)T,L{k), yielding after some 
calculations 

=l9'T'(^N^ + C.Y.^{^)^ xigArn,}), (15) 

00 

with the same value for the transverse and longitudinal polarizations, giving rise to a 
mass gap, i.e. uj{k — > 0) > 0. In general, x{9: {^g}) depends on a vector {mg} with 
components mg. The function x((7,m) = x((yf, {m, m}), depicted in the left panel of 
figure [3]for one independent mass m = rUg and Nf = 2, is implicitly defined by (|T5l) with 
the numerically determined value of ujpi = uj{k —>■ 0) on the l.h.s. It approaches unity 
for small values of g in agreement with the HTL result. At large values of the coupling 
g, x(5f,m) rapidly drops, making the plasma frequency soft. In addition, x((yf,m) is 
almost independent of m in the range of interest. From this follows that the mass 
dependence of the plasma frequency is essentially contained in the function J{m/T), 
which is depicted in figure [3], right panel. For large values of m/T, J{m/T) behaves 
like exp(— m/T) indicating the decoupling of the heavy quark sector from the thermal 
bath, while for decreasing values of m/T, JimjT^ approaches unity with m/T — > 0. 
Thus, in the limit (7 — > and m — or T ^ 00, equation (fT5i) gives the usual HTL 
gluon plasma frequency CJ^i = ^g'^T^{Nc + C2Nf). 

In the asymptotic region, k 00, the transverse gluon dispersion relation is given 
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Figure 2. Transverse (left) and longitudinal (right) gluon dispersion relations 
represented as Qt,l = {^t,l ~ ^^)/'^g functions of scaled momentum and scaled 
quark mass for g = 0.3, 1.0, 3.0 from top to bottom. For a better presentation, the 
transverse and longitudinal quantities are considered from different perspectives. Faint 
isolines are given at Gt.l = 0.1, 0.2, . . . , 0.9. 



by the gauge invariant result u)^ = k"^ + with the asymptotic gluon mass 



Re UUk, k) = \g'T' |^iV, + ^2 

oo 

e 1 

v/rTei+exp(xv/rTe)' 
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9 m/T 

Figure 3. Left panel: Function defined in (fTSll . showing the non-trivial 

dependence of the gluon plasma frequency as a function of the coupling g. On the 
displayed scale, curves for < m/T < 2 are on top of each other. Right panel: 
Functions J (dashed curve) and I (solid curve) as defined in (jl6p and (jlSp , respectively. 



The function 1 is depicted in the right panel of figure [31 The power expansion for small 
values oi m/T (cf. Appendix A) reads 

OO 

X [j) = l+«2 [j) +«L [j) In (-) +a, (-) (-) (19) 

with coefficients a„ listed in table IAT] relegated to Appendix A. This expansion of X is 
convergent in the region m/T < vr. Note the term oc ai resembling chiral logarithms. 
Equations f|T7j) and f|T8|) highlight an important feature occurring in the gluon sector, 
here, for the example of asymptotic dispersion relations. As the quark loop contribution 
is added to the gluon and ghost loops, the mass dependence is not so striking in the range 
of m we are interested in. For instance, for Nj = 2 the term in brackets in (fTTj) varies in 
the range between 4 and 3 while going from m = to m ^ oo, evidencing the dominance 
of the strong non-Abelian self-coupling compared to the quark loop contributions in the 
thermal masses. In a QED plasma there is, of course, a much stronger dependence of 
the thermal photon mass on the electron mass m^. 

The heavy-quark expansion for m/T ^ 1 yields, as evident from (1A.2I) in Appendix 

A, 

_ /m\ 12 Inm m , , 

As I{m/T) becomes exponentially small for large values of m/T, the heavy quarks 
decouple also in the asymptotic momentum region from the thermal bath (i.e., they are 
exponentially suppressed) . 
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3.1. Quark self-energy 



The one-loop quark self-energy represented in the imaginary time formalism reads 



(21) 



(27r)3^^ 

with loop momenta P and Q — K — P. The general structure of the self-energy reads 

T, ^ afC + byt - 7oK^> k) - 7ka(a;, k) - c(u;, k) (22) 

with the medium 4- velocity and b — u!a-\-b. After performing the Matsubara sum and 
continuing into Minkowski space-time one finds for the three independent self-energy 
functions 



a{u!, k) 



27r^k'^ 



dp 



I . ( ^ ^ hp , + dp , 
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with Cp = (N^ — l)/{2Nc), and the abbreviations read 



ap = 



Ob 
hp 



k'^ — vn? — uj"^ ± 2€pUJ — 2pk 
^2 _ _ _|_ 2epUj -\- 2pk ' 
k'^ -\- — u;"^ ± 2puj — 2pk 

_|_ _ (^2 _|_ -)- 2pk ' 

A;^ — — o;^, c?f = SepO;, 
k"^ -\- n? — up' ., ds = 2pu . 



(23) 
(24) 

(25) 

(26) 

(27) 

(28) 
(29) 



3.2. Quark dispersion relations 

The dispersion relations are determined by the poles of the resummed quark propagator 
<S, given by the Dyson equation — S'^^ -\- E, leading to 

= 7o(cj + 6) -7k(l + a) - (m + c). (30) 

With the definitions of projectors relegated to Appendix B, this can be rewritten yielding 

(31) 



r(rza; — 8) r{nuj + 8) 
The quark (plasmino) dispersion relation is obtained as the positive energy solution of 
nUq — 8 = Q {nujp -\- 8 = 0). The dispersion relations for the quasi-particles may be 
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parametrized as cu^ = fc^ + + m'jj^g and = k'^ + m'jj^p for quarks and plasminos, 
respectively. The functions JF^ ^, which non-trivially depend on k and m, encode the 
effects on the excitations induced by the thermal medium with temperature T. For fixed 
k and m these functions multiplied by m'j serve as effective thermal mass parameters. 
J-'q = accounts for a simple on-shell dispersion relation, while J-'p = represents the 
light cone. The quantity m'j = ^Cpg'^T'^ is the usual HTL quark plasma frequency, 
setting a typical energy scale in these studies . 




Figure 4. Quark (left) and plasmino (right) dispersion relations represented as 
J^q — {uj^ — fc^ — m'^)/m'j for regular quark excitations and J^p = {cOp — k"^) / m'j 
for plasmino excitations as functions of scaled momentum and scaled quark mass for 
g — 0.3, 1.0, 3.0 from top to bottom. For a better representation, the regular quark 
and the plasmino quantities are considered from different perspectives. Faint isolines 
are for Tq^p = 0.1, 0.2, . . . 1.9. The representation of J-p in the right panel ends where 
numerically no real solutions of the dispersion relations can be found. 
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A survey of this representation of the dispersion relations is exhibited in figure HJ 
where the temperature dependent parts JF^ and Tp for regular quark and collective 
plasmino excitations are shown in the left and right columns, respectively. Again, 
different cuts allowing for a better quantitative representation are relegated to Appendix 
C (cf. figures [C3] and [C4l) . Further examinations of the one- loop self-energies show that 
these exhibit non-zero imaginary parts, both, below and above the light cone which are 
particularly large below the light cone. 

Considering, instead, ujp = k'^ + m'jj-'p as a function of k"^, the plasmino branch 
exhibits a local minimum which, however, vanishes for m > gT. This local minimum 
gives rise to van Hove singularities in certain emission rates |38j . 

In addition, for any value of m > 0, the collective branch completely disappears 
from the spectrum for momenta larger than a certain critical one. Furthermore, 
the plasmino branch lies energetically below the regular quark excitation branch 
(Surprisingly, the lattice QCD calculations [25] show an opposite order at k = 0). 

As mentioned above, with some replacements (here, one has to use Cp —>■ 1 in 
addition) the electron excitations in an electron-positron QED plasma are also described. 
Particularly interesting would be an experimental verification of the purely collective 
plasmino excitations. 



3.3. Quark plasma frequency and asymptotic behaviour 

For vanishing momenta k ^ 0, the plasma frequencies for quarks (+) and plasminos 
(— ) uj± = uj{k —>■ 0) are implicitly given by 

uj^^m+ (bTc) =0, (32) 

where we used {ka)\f^^Q = due to rotational invariance. The long- wavelength limit 
of the quark dispersion relation is therefore independent of the function a. An explicit 
calculation of the various limits yields the implicit equation 

oo 

n ^ , d'^^F f dpp \ pep 

= uj^^rn + ^ J ^ n^(ep)-^(^±Tm) 



nsip) 



+ ^^[{p'Tmz±) + z±{uj±Tm)) (33) 

with z± = (m^ —uj^)/2u}±. The splitting of the two branches is of order m; in the chiral 
limit, m ^ 0, both branches meet at /c = 0. In the limit (7 -C 1 and m < gT, one 
recovers a formula given by Pisarski [30] 

a;± = - i^-^m? + 4mj ± , (34) 

and the difference of the plasma frequencies of quarks and plasminos becomes a;+ — uj^ = 
m. Thus, the regular quark has a higher excitation energy than the plasmino. In the 
limit m = the high temperature result uj± = mj is obtained. 

In the asymptotic region, the momentum k is the largest scale, allowing for some 
approximations to obtain an explicit analytic expression for the asymptotic dispersion 
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relation of the regular quark excitation. In this case, the equation nu — S = can be 
transformed to yield 



uj"^ = k'^ + + 2 [k'^a + kb 



, (35) 

fe— >oo 

which is completely independent of the self-energy function c. The asymptotic quark 
dispersion relation can then be written as 



e + m' + 2Ml, = ImJ (j (^) + 2) (36) 

with the integral X given in (ITSl) and the same expansion as in (fT9l) . For m ^ 0, one 
finds — i> nij. 



4. Gauge dependence 

As the HTL approximation has been proved to be gauge independent [33], one may 
look for regions in parameter space where the self-energies in one-loop approximation 
coincide with HTL results. The real parts of the HTL and one-loop self-energies coincide 
for small values of uj and k, of course, but for larger values of u and/or k both may 
deviate noticeably in general. However, near the light cone, i.e. for uj ^ k, they are 
approximately equal , and on the light cone, u = k, HTL and one- loop self-energies 
coincide, cf. [HI [39]. The self-energies slightly above the light-cone determine the 
excitation energies of the quasi-particles in a wide parameter range and are, therefore, 
relevant for the quasi-particle model [3, [9] . 

The gauge independence of the asymptotic gluon mass can be shown by noting 
that the transverse parts of the Yang-Mills contributions to the gluon self-energy at the 
light cone Il^^{k,k) are exactly the same in the one-loop and HTL approximations. 
Furthermore, the quark contribution H^ is inherently gauge independent at one-loop 
level. Hence, the asymptotic thermal gluon mass is a gauge invariant quantity. 

To check further for the gauge dependence of the above one-loop results in Feynman 
gauge - here in particular we focus on the asymptotic mass of regular quarks - we 
compare them with explicit calculations in Coulomb gauge. For large momenta, the 
regular quark excitation energies in Feynman gauge (u) and Coulomb gauge (uc) are 
related via 

ujQ{uj,k) = u'^^u, k) + Ac{uj, k), (37) 
00 

A ( M s'Cf r ( P + p^ k^ + p^ + 2kp \ 



with lim Ac{uJ, k) — * 0, thus proving the coincidence of the dispersion relations in both 
gauges in the asymptotic region. 



5. Summary 



In summary, we survey the quark mass dependence of thermal one-loop self-energies 
and emerging dispersion relations of quarks and gluons in QCD in Feynman gauge. The 
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motivation and the related focus of our presentation are given by upcoming lattice QCD 
results on spectral properties of quarks and gluons in a hot and deconfined medium and 
the need of chiral extrapolations in bulk properties of the quark-gluon plasma. While 
the results of ^25] are at variance with the energetic ordering of quark and plasmino 
excitations from perturbative QCD, they otherwise confirm the existence of a mass gap of 
quarks in line with basic assumptions of the successful phenomenological quasi-particle 
model [3 H] . The present analysis challenges the chiral extrapolation of the equation 
of state performed in [10], as the important term oc mM_|_ in the quark excitation 
dispersion relation is not supported by our one-loop results. This may be considered 
as further hint to strong non-perturbative effects in the quark-gluon plasma in a range 
accessible in present and future heavy-ion collisions. Indeed, the asymptotic Dyson- 
Schwinger type approach outlined in Appendix D, points to severe corrections to the 
one-loop results in the strong coupling regime for fermion masses of the order of the 
temperature. Progressing lattice QCD results are required to reveal the fundamental 
excitations of deconfined matter. 
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Appendix A. Expansion of the function X 

To derive the expansion in ( fT9i) of the function X given in (fT8|) we rewrite the expression 



using (cp + 1) 



1 



oo 





and obtain 




(A.2) 



with X = m/T. For small quark masses, x 0, one may represent I as 




(A.3) 



with the Mellin transform 



oo 




(A.4) 
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of the first modified Bessel function of the second kind Ki. 

Exploiting = (1 ^ 2~'*) C^{s + 1), with the Riemann zeta function C, 

we find as intermediate step 

c+ioo 



X(x) 



12 1 

TT^ 2'Ki 



ds 



X 



l-s 



1 



r 



s + l 



which can be evaluated with a suitable contour and appropriate c 
J = 5^Res/(.)|_^, 



2 

to yield 



f{s) 



6 X 



2 Ol-s 



7r2 2 



1 



C(5 + i)r 



1 



r 



1 



, (A.5) 



(A.6) 



(A.7) 



Further evaluation leads to the series in (fT9|) with coefficients given in table lATl 
Although the numerical value of the coefficients rapidly approaches zero, this expansion 
has a finite radius of convergence, which is 



lim 

n— >oo 



an+2 



TT 



n > 4. 



Table Al. The coefficients of the chiral expansion in (fT9|) . is the Euler-Mascheroni 
number and C stands for Ricmann's zeta function. 



coefficient analytic expression numerical value 



a2 


-;^(ln7r + i-7B) 


« -0.32449 


CtL 


3 


w 0.15198 


a4 


-Ti-CO) 


« -1.619- 10-2 


ae 


1287r'iC(5) 


« 7.836 • 10-4 


as 


1905 


-4.943- 10"^ 



for n = 4, 6, 8, . . . 



Appendix B. Quark propagator 

The explicit breaking of chiral symmetry drastically modifies the structure of the quark 
propagator compared to the usually treated case, where all current quark masses are 
set to zero. As a consequence, the usual projectors, which allow for a distinction of the 
different quasi-particle states, are not applicable. Thus, one has to examine in some 
detail the spinor structure of the propagator, which provides projectors to the relevant 
subspaces of the physical (quasi-particle) excitations. For the self-energy we employ the 
decomposition given in (1221) . Utilizing Dyson's equation yields the general structure for 
the inverse quark propagator 

S-^ = r {ncu-fo - jk - M) (B.l) 
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with the three self-energy functions r = (1 + a), n = l + h/uor and }A = (m + c)/(l + a), 
depending on uj, |k| and m. Although there is an explicit breaking of chiral symmetry 
due to the finite current quark masses, chiral symmetry is assumed not to be broken 
spontaneously at high temperature, thus, lim A^(ct;, fc, m) = 0. In contrast to the 
vacuum, an additional function n appears, which may be interpreted as an analogue 
to the refraction index in optics [H]. The appearance of this function is due to the 
breaking of Lorentz invariance by the term oc ^ in (1221) . 

The spinors ip representing the quasi-particle states are supposed to be solutions of 
the Dirac equation 

(ncJ7o - 7k - M)i)ik) = 0. (B.2) 



Non-trivial solutions are found only for n^cj^ = S"^ with S = -\-y/k^+Ai^. This is an 
implicit equation for the excitation energies uj of the system, leading to u = ±S/n. A 
Hamiltonian may be defined by 

n = -i^olk + ^oM) (B.3) 
n 

allowing for rewriting (IB. 21) as Tiipik) = ujipik). Let Ua be a solution of this 
equation, then we have UaUJa = ±Sa = ±a/A;^ + Ail, where Ua = n{uja,k,m) and 
Aia = -Mi^^a, k,^), and the states can be classified by the eigenvalues of the sign 
operator A^^ = uTi/S, which obviously commutes with Ti and has eigenvalues At^ = 
sgn(nci;) = ±1. 

In the limit m — 0, the operator A^^ becomes the operator of chirality times helicity, 
A = 707k, where k = k/|k|. Thus, the states with eigenvalue +1 have the same value 
for chirality and helicity, whereas the states with eigenvalue —1 have the opposite value 
for chirality and helicity. This simple interpretation fails for m 7^ 0, as chirality is not 
a conserved quantum number in this case. 

The projectors that decompose the Dirac structures according to the physical 
(quasi-particle or collective excitation) states are induced by the operator A-j^ and read 



^ , 7k + A^ 

1 ±7o- 



(B.4) 



satisfying Vtm'^tm = and V^V^,^ = 0. Then, %^ projects on all 

possible states with sgn(?T,u;) > (sgn(nci;) < 0). In the chiral limit m — > 0, they take 
the form = |(1 ± 707k), and for A; — > one obtains = |(1 ± 70), which are 
often employed decompositions. Our more generally decomposed quark-propagator for 
m > takes eventually the form of (1311) . 



Appendix C. Cuts through the dispersion relations 

We present here various cuts of m = const, g = const and k = const to quantify 
the dispersion relations and expose the relevant dependencies. Figure ICll exhibits Qt,l 
from figure [2] in section [2lB as a function oim/T for various fixed values of momentum 
k (different curves) and couplings g (different panels). The mild mass dependence 
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Figure CI. Cuts through the dispersion relations exhibited in figure [2] as a function 
of m/T for k/gT — (red sohd curves), 1 (green long-dashed curves), 10 (blue 
short-dashed curves) and the asymptotic region (fc ^ oo, black dotted curves) for 
g = 0.3, 1.0, 3.0 from top to bottom. Left: Transverse gluon branch Qt, right: 
Longitudinal gluon branch Q^. This figure exposes the dependence on m. 



mentioned in section OB is clearly seen in this representation. For longitudinal modes, 
uj{k) approaches the light cone for larger momenta. In addition, at small k the 
dependence on the coupling g is strong for all considered values of m, while at large 
momenta the coupling strength is of minor importance (cf. transverse modes). Note the 
kinky structure for small values of k. The momentum dependence of Gt,l is displayed 
in figure \C2\ for various fixed values of m/T (different curves) and various fixed values 
of the coupling g (different panels). Again, the weak m dependence is clearly visible. 
Note also the fairly strong dependence on the coupling strength. 

Figures IC3I and IC4I exhibit various cuts through the dispersion relations of quark 
excitations shown in figure H] in section [31 B in an analog manner as presented for 
the gluons above. Figure IC3I exposes the strong quark mass dependence for various 
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Figure C2. Cuts through the dispersion relation exhibited in figured] for m/T — 
(red solid curves), 0.1 (green long-dashed curves) and 1 (blue short-dashed curves) 
and for g = 0.3, 1.0, 3.0 from top to bottom. Left: Transverse gluon branch, right: 
Longitudinal gluon branch. This figure exposes the momentum dependence. 



fixed values of momentum (different curves) and coupling strengtlis (different panels). 
Interesting is the non-monotonic beliaviour at small momenta which is most pronounced 
for weak coupling. Striking is the disappearance of the plasmino branch at large 
momenta but also for smaller momenta and increasing mass, (see right panels). The 
plasmino branch (if there is any) persists for larger values of m the larger the coupling 
g is. The energy splitting with changing coupling is most severe for small momenta. 

The momentum dependence of jFg and JF^ is exhibited in figure IC4I for various 
fixed values of m/T (different curves) and various values of the coupling g (different 
panels). For normal quark excitations Tq rises monotonically with increasing momenta 
(left panels), while for plasmino excitations J^p drops with increasing momenta (right 
panels). For strong coupling g, the plasmino excitations become independent of the 
quark mass (right bottom panel), while in the weak coupling regime a severe mass 
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Figure C3. Cuts through the dispersion relations in figure S] as a function of the 
scaled quark mass for k/nif = (red soHd curves), 1 (green long-dashed curves), 10 
(blue short-dashed curves) and asymptotically large values (M^ , black dotted curves) 
and for g = 0.3, 1.0, 3.0 from top to bottom. Left: Normal quark branch (J\j), right: 
Plasmino branch (J-p). This figure exposes the dependence on m. 



dependence is visible (right top panel). Note also the approximate mass independence 
of J-'q for small quark masses in the strong coupling regime (left bottom panel). 

In the quasi-particle description of QCD thermodynamics (cf. [3 El), equation f|T7j) 
is used with ~^ i-^- the quark mass dependence in the asymptotic gluon 

dispersion relation is disregarded. Inspection of the left panels of figure ICll supports 
this approximation: The energy of transverse gluon excitations (say, with momenta 
k ~ gT) increases only by a tiny amount, in particular for larger g, when decreasing 
the quark mass from m/T = 0.4 (as used in [31] for "light quarks") to the chiral limit 
m = 0. For larger values of k this increase is still less than 4 %. Curiously, in the chiral 
limit, the one- loop transverse gluon excitations obtain a larger thermal mass, thus, the 
according entropy carried by these modes is reduced. 
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k/lTLf k/lTLf 

Figure C4. Dispersion relations Tq^p vs. scaled k as in figure 2] but for various cuts 
with m/T — (red solid curves), 0.1 (green long-dashed curves) and 1 (blue short- 
dashed curves) for g = 0.3, 1.0,3.0 from top to bottom. Left: Normal quark branch, 
right: Plasmino branch. This figure exposes the momentum dependence. 

Concerning the asymptotic (large momenta) regular quark dispersion relation, a 
bi-linear term 2mM+ is employed additionally in the quasi-particle model [7[ [9] which is 
absent in ( !36l) . (In fact, also any quark mass dependence in is disregarded in [71 [9], 
thus, reducing to m^.) The term 2mM+ used in [71[9] was inspired by [30], mimicking 
the quark mass dependence of regular quark excitations at small momenta and small m 
and g (see also red curves in the left panels of figure IC3I for small m and g) . There, in 
particular for small and moderate couplings, drops by a significant amount (about 
50 %) when decreasing m from m/T = 0.4 to the chiral limit m = 0. As a consequence, 
the entropy attributed to such excitations increases when lowering the quark mass. In 
[40] (last reference) it was shown that with such an ansatz for the quasi-quark dispersion 
relation by adjusting the quasi-particle model parameters to lattice QCD results [3T] 
(first reference) using m/T = 0.4 for "light quarks" good agreement with lattice QCD 
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thermodynamics [12] for " almost physical quark masses" was found when extrapolating 
to the corresponding quark mass values. 

In contrast, based on the one- loop considerations presented in this work, the thermal 
mass contribution Tq of regular quark excitations increases with decreasing m for large 
momenta [Tq = 2M^/m'j resembling the behaviour of the function I with m), as seen 
in the left panels of figure IC3I If such a quark mass dependence (i.e. equation fl36l) ) 
would have been implemented in the quasi-particle model, the quark mass extrapolation 
reported in [lO] would fail to reproduce the corresponding lattice QCD results [12]. This 
highlights the subtle role of the actual dispersion relation employed in the quasi-particle 
model, as mentioned already in [13] and points to the necessity of a term oc mM+ 
of sufficient strength in order to account for the quark mass dependence in the non- 
perturbative regime. 

It happens that for momenta k/T ~ C>{1) and small m, the small-momentum 
ansatz amM+, where a ^ 2, represents a better approximation to the mass dependence 
of the regular quark dispersion relation than the asymptotic form in (l36l) (see left panels 
of figure [C3| . This motivates to some extent the dispersion relations employed in the 
model and the performed chiral extrapolation in [IQ]. (Note that the quasi-particle 
model [3 [9] does not involve plasmino and plasmon excitations; their contributions to 
the entropy density, for instance, are found to be numerically small [44J.) 



Appendix D. Asymptotic Dyson-Schwinger approach in Abelian gauge 
theory 

In order to estimate higher-loop order correction effects on our one-loop results, we 
contrast the one-loop results with corresponding calculations obtained in a Dyson- 
Schwinger type approach. In the following, we focus on the fermion mass dependence 
in the asymptotic thermal mass expressions. The Dyson-Schwinger equations (DSE) for 
the 2-point functions, here presented for an Abelian gauge theory, read 




for fermions and 




iv-')^u = iv,')^u + n^. (D.2) 

for bosons, where the self-energies S and 11^,^ are functionals of the dressed propagators 
S and V^^ (indicated by fat blobs) as well as of the dressed vertex-function (also fat 
blobs). Thus, these expressions couple to higher order DSE for higher n-point functions. 
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Figure Dl. The mass dependence of the scaled asymptotic thermal masses and 
vriF for Abelian gauge bosons (left panel) and fermions (right panel), respectively. The 
self-consistent solutions of the gap-equations (|D.5P and (|D.6|) are represented by solid 
and dashed lines for various values of the coupling g. The black dotted lines show the 
one-loop results given in (fT7|) and ([551) . respectively, using the substitutions described 
in the text to go from non-Abclian (QCD) to Abelian (QED) gauge theory. In the 
notation of the preceding sections this would be m^/m^ for the gauge boson and 
M^/(2m^) for the fcrmion. 



In order to calculate the asymptotic masses from the DSE, we use momentum 
independent asymptotic thermal mass expressions as ansatz for dressing the propagators, 
because the asymptotic thermal masses result from mutual scatterings among hard 
plasma particles with no coupling to the soft momentum excitations. (This ansatz is in 
an analog spirit as the approach in [15].) Then, the propagators read in Feynman gauge 

fi — mp 



where we only dressed the transverse part of the gauge boson propagator [16]. The 
Ward-identity, dS^^ jdK^ = ■y'^ suggests to use the bare vertex, giving rise to a natural 
decoupling from the higher order DSE. We focus, here, on the poles of these propagators 
at hard momenta by solving the DSE in ( ID. II) and ( ID. 21) in the asymptotic region. As 
a result, by using only the temperature dependent parts of the self-energies, we find for 
the gap-equations 



m| = ^iFimp) (D.6) 





2 
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with 



oo 




dk — np^ep), 



(D.7) 







oo 




dk UBitB), 



(D.8) 







where e\ = k"^ + m^, e\ = k"^ + rn^p. Self-consistent solutions of these equations for m\ 
and m\ are presented in figure [Dll For small values of the coupling, say g = 0.03 (cyan 
long-dashed curves), the self-consistent solutions confirm the mass dependence in the 
one-loop results quite well. For larger values of the coupling, some deviations between 
the two approaches emanate, in particular in the fermionic sector, where already for 
g > 0.3 the deviation becomes visible. In the strong coupling regime {g > 1) the 
difference becomes 50% and larger. The coupling dependence of the self-consistent 
fermionic mass at m = shows a similar behaviour as was found in [23] insofar, as in 
both approaches the self-consistent mass is much smaller than the corresponding HTL 
or one-loop values, although [23] calculated the long-wavelength limit of the fermion 
dispersion relation while we concentrated on the asymptotic limit. 
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